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g ; Abstract 
o : 

CN , We obtain tree-level estimates of various differential branching ratios of heavy baryon decays 

I with massive sterile neutrinos in the final state. Heavy sterile neutrinos can be searched for 
Q ' in heavy baryon decays with missing mass as a signature as well as in fixed target experiments, 
^ ■ where heavy baryon decays contribute to sterile neutrino production, with sterile neutrino 
decays as a signature. Generally, charmed baryons are found to be less promising than charmed 
r-| ! mesons, in contrast to b-hadrons. In the latter case, branching ratios of beauty mesons and 
baryons into sterile neutrinos are of the same order. As a consequence, at high energies beauty 
Oh| baryons give contribution to sterile neutrino production comparable to the contribution of 
beauty mesons (up to about 15%). Experimental limits on active-to-sterile mixing are quite 
strong for neutrinos lighter than D-meson but for heavier neutrinos they are weaker. As an 
y—i , example, for neutrino masses in the range 2 GeV < m^^ ^2.5 GeV, current data imply that 
^ [ the bounds on A^-hyperon branching ratios into sterile neutrinos are Br(Afe ^ Ac + e~ + Vx) ^ 
§ ■ 1.3 ■ 10-5 - 1.7 ■ 10-6 and Br(Ab A^ + /x- + p^) < 3.9 ■ 10-^ - 1.4 ■ 10-^. 

\o : 

p ; 

O 1 Introduction 

oo ! 

Q , Sterile neutrinos are introduced in particle physics to explain neutrino oscillations: they provide 
j>! I the active neutrino sector with masses and mixing. The inferred mass and mixing patterns in 
the active neutrino sector are severely constrained by neutrino experiments whereas resulting 
limits on parameters of the sterile neutrino sector are strongly model-dependent. In particular, 
sterile neutrino masses can be in the GeV range and this fact does not necessarily imply very 
small values of mixing angles [T]. In this case sterile neutrinos can be searched for in particle 
physics experiments, and here we discuss in detail sterile neutrino production in baryon decays. 

For this study we assume that some sterile neutrinos have masses in the GeV range, and 
consider active-to-sterile mixing angles as free parameters constrained from direct searches for 
sterile neutrinos. This model-independent approach yields most conservative constraints on 
the sterile neutrino production rates. While chosen ad hoc, the GeV range for sterile neutrino 
masses can be preferred in some models, with z^MSM [21 E] serving as an example. 

The richest source of sterile neutrinos of masses in the range we consider are weak decays 
of heavy mesons created in beam-beam and beam-target collisions. Recently, these decays 
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were revised in Ref. If sufficiently light, sterile neutrinos are produced mostly in tt- and 
i^-meson decays. For heavier neutrinos, leptonic and semileptonic decays of charmed mesons 
are the most relevant sources of sterile neutrino production. Finally, in models with neutrinos 
heavier than D-mesons but still in the GeV range, decays of B-mesons dominate. In all cases 
additional contribution comes from decays of baryons. For relatively heavy neutrinos the largest 
contribution from this baryon channel is due to semileptonic decays of charmed and beauty 
baryons, among which the most promising are A^- and A^-hyperons because of the largest 
statistics collected. 

In this paper we present explicit formulas for differential branching ratios of baryon semilep- 
tonic decays with massive sterile neutrinos in the final state. These formulas can be used to 
estimate the contribution of the baryon channel to the spectrum of sterile neutrinos produced 
in fix-target experiments, where sterile neutrino decays can be searched for. Likewise, the ob- 
tained formulas are relevant for searches of heavy baryon decays into sterile neutrinos which 
can be performed at operating and future B-factories, Tevatron and LHC 

2 Probabilities of baryon decays 

We start with the semileptonic decay of Ac-hyperon. Ac-hyperon weakly decays into A-hyperon, 
positron and electron neutrino (as well as positive muon and muon neutrino). Let sterile neu- 
trino Kiix with electron neutrino or niuon neutrino and Oi^ {I = e, ^i) is the corresponding 
mixing angle. Then the amplitude of the process Ac ^ A + + z^^ can be written as follows, 

G 

M = -^Vcs sin 



where P and (P — q) are 4-momenta of A- and Ac-hyperons, Ma and Ma, are their masses, 
G is the Fermi constant, Vcs is the element of the Cabibbo-Kobayashi-Maskawa matrix. The 
dimensionless form factors /i, /2, /s, gi, g2, gs entering Eq. ([T]) parametrize the matrix element 
of the relevant hadronic current = 07,^(1 — 75)^ between real A- and Ac-hyperons. These 
form factors are functions of q^. 

The differential decay rate is given by 
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P, ky^, ki denote the 3- momenta of A-hyperon, sterile neutrino and charged lepton, respectively, 
i?A, Ey^, El are their energies, is the squared amplitude, averaged over spins of the initial 
baryon and summed over spins of final particles. 

The direct calculation gives the following expression for |A1| : 

\M\^ = 8G^\Vcs\'' sin^ Oi, \if^ + gl){APk,^ Ph - 2Pk,^ qh - 2Pki qk,J 



-2{f^-gl)Mj,Mj,X^k 



^(/i/2 + 9ig2)ik,MPq - + qKAPki - gk) + gkiPk,^ - qk,j) 

2(/i/2 - 9i92)iPki qk^^ + Pk^^ qh + Pq k^ki) 
^^^^(4Pfc,, Ph - APq{Ph qk,^ + Pk,^ qh) 

- 2k,M{'^{Pq? - q'Ml - Pq q') + 2qk,^ qh [Pq + Ml) 

- 3Mlq^ k,^h + A{Pqfk,^h - Pq q^ k,^h) 

if 2 - 9l){q^ KM + '^qK^ qh) 



+ ^(/i/3 + gigz)M^M^{2qk,SPh - qh) + 2qh{Pk,^ - qk,^) - 2k,MPq ' q')) 
+ (/1/3 - 9i93)i2Ph qk,^ + 2Pk,^ qh - 2Pq k,^h) 
+ 2l^I^^q^{Phqk.^ + Pk.^qh) 

(2qk,^ qh - q^ k,ki) {Pq - M^{M^^ + Ma))) 



Ml 



{2qk,^ qh - q^ k,^h) {Pq + Ma(Ma, - Ma))) 



Ml 



+ A{Ph qku^-Pk^^qh) 



(3) 



Integrating Eq. ([2]) over 3-momenta of sterile neutrino and charged lepton and over all 
possible directions of the outgoing baryon, one gets for the differential decay rate 



dT _G'\V,s\' sin2 ^El - Ml q' - 2q^{ml + mf) + (mf - 



dEA 64773 q4 y q4 

^ {2q' - q\ml + m^) - K - ((/i + g^mW - l2Pq q' + 8(Fg)2) 
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- 12{f^ - gl)M^M^£ + 24^(/i/2 + gig2){Pq - q')q' - 2A{hh - g,g2)Pq q' 

- ^^{12Pq q' - 8{Pq)V + AMlq') - Afl ' 9l)q') 
+ ((m^ - ml^f - q'iml + m^)) {{f^ + gl){AMlq^ + APq q' - 8(Fg)2) 

- 4(A^ - gl)q'Mj,M^^ - 8^(/i/3 + g,g3)Pq q' - 8(/i/3 - gm)Pq q' 

+ - ^a(Ma^ + Ma)) + ^^q\Pq + Ma(Ma, - Ma))) 

Ac Ac 

where Pq and q^ are functions of Ep^: 

Pq = Ml- EaMa^ 

and 

q^ = Ml + Ml - 2M^E^. 

Semileptonic decays of any other baryons {Bi 32 + 1 + ^) are described by similar formulas 
with obvious replacements: Ma^ — > M^^, Ma M^j, |V^sP — ^ l^iP) where Vij is the relevant 
element of the Cabibbo-Kobayashi-Maskawa matrix. 

By integrating Eq. (jll) over the final baryon energy one obtains the decay rate F as the 
function of sterile neutrino mass m,^^. The corresponding branching ratios of the decays — >■ 
A + + z/j; and A^ — > Ac + /~ + i^i; are presented in Fig. [H charged lepton is considered to be 
massless (electron or positron). Results of numerical calculations for the differential branching 
ratios, which give differential spectrum of outgoing baryons, are presented in Fig. O 

In all numerical calculations we use the form factors in the dipole approximation [5]: 

fM') = -7-^^, (5) 

1-4 

9i{q ) = — 2> (6) 

1-4 

where for charmed baryons mv = 2.11 GeV, tua = 2.54 GeV, and for beauty baryons my = 6.34 
GeV, niA = 6.73 GeV. Values of /j(0) and (?i(0) for the transitions A^ Ac and Ac — A are 
taken from Ref. |5]. They are summarized in Table [H Values of /2(0), /3(0), g'2(0), gsi^O) have 
the opposite sign as compared to Ref. [5] because of the different parametrization of matrix 
element ([T]). 
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Form factors 


Afc Ac 


Ac ^ A 


/i(0) 


0.53 


0.29 


/2(0) 


0.12 


0.14 


/3(0) 


0.02 


0.03 


^7i(0) 


0.58 


0.38 


^?2(0) 


0.02 


0.03 


^?3(0) 


0.13 


0.19 



Table 1: Form factors for the transitions Ac — > A and A;, Ac adopted from Ref. |5] and used 
in numerical calculations. 
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Figure 1: Branching ratios 
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for baryon decays: a) Ac ^ A + + z/^.; b) A5 — > Ac + / \Vx- 



The differential rate describing the differential spectrum of outgoing charged leptons, 



is obtained by integrating Eq. 
dV 1 



over P, ky^. The final expression is: 



dE, 647r3MA 
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where 



p^ = Ml^ + 7ii?i-2M^Eu 



9 is the angle between charged lepton and outgoing baryon in the center-of-mass frame of sterile 
neutrino and outgoing baryon. The results of numerical calculation are presented in Fig. [3l 
Differential spectrum (JTj) can be used in searches of heavy baryon decays with sterile neutrinos 
in the final state. 

The differential rates are given by the same expression with obvious replacement: 
l] is now the angle between neutrino and outgoing baryon in the center-of-mass frame 
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Figure 2: Differential branching ratios '^^d(E^^G<^f functions of final baryon energy for baryon 
decays: a) Ac ^ A + /"*" + z/^, various plots correspond to various sterile neutrino masses: from 
top to bottom my^=Q, 0.2, 0.4, 0.6, 0.8, 1.0 GeV; b) Kb ^ k^ + l^ + t^x] rn^,=Q, 0.5, 1.0, 1.5, 
2.0, 2.5 GeV. 
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Figure 3: Differential branching ratios '^^fiEiJccY)^ functions of charged lepton energy for 
baryon decays: a) Ac — > A + /"*" + z/^.; b) Af, — Ac + /~ + Ux- Sterile neutrino masses are the 
same as in Fig. [21 Larger branching ratios correspond to smaller neutrino masses. 



of lepton and outgoing baryon. The numerical results are presented in Fig. H] and are relevant 
for searches of sterile neutrinos in baryon decays, since describes the spectrum of produced 
neutrinos. 
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Figure 4: Differential branching ratios '^^^^'^/gcv')'' functions of sterile neutrino energy for 
baryon decays: a) Ac — A + + z/^.; b) A^ Ac + /~ + h'x- Sterile neutrino masses are the 
same as in Fig. [21 Larger branching ratios correspond to smaller neutrino masses. 



3 Discussion 

Our results for the baryon branching ratios into sterile neutrinos can be compared to similar 
results for heavy meson decays, presented in Ref. In the case of charmed hadrons the general 
conclusion is that baryon branching ratios are always significantly lower than branching ratios 
of charmed mesons. As an example, for niK ^ rriiy^ < 1 GeV, where rriK is K-meson mass, 
one estimates ^-j — 10. Branching ratios of leptonic meson decays are even larger. 

Hereafter, meson branching ratios are calculated with the use of formulas presented in Ref. [3] 
and form factors presented in Ref. [6]. Hence, charmed baryons are less promising for sterile 
neutrino searches in comparison with charmed mesons, since larger statistics is required to reach 
the same level of statistical sensitivity to active-to-sterile mixing angles. From the analyses of 
baryon and meson decay modes, presented in this work and in Ref. [1], respectively, and with 
account of baryon fraction in c-quark hadronization [7j one can estimate the contribution of 
baryon channel to sterile neutrino production in proton beam-beam and beam-target collisions 
at the level of 2%. 

In beauty sector the situation is quite the opposite: baryon branching ratios are somewhat 
larger than branching ratios of mesons. As an example, for 2 GeV < < 3 GeV we obtain 
bJ{a^^a ^ ) — 0.5. Note that B ^ D* + 1 + dominates the sterile neutrino production in 
meson channel. Thus, one expects that searches for sterile neutrinos in baryon decays should 
be competitive with similar searches in meson decays. From the analyses of baryon and meson 
decay modes, presented in this work and in Ref. |3] , respectively, and with account of baryon 
fraction in b-quark hadronization [8j one can estimate the contribution of baryon channel to 
sterile neutrino production in proton beam-beam and beam-target collisions at the level of 15%. 
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The absolute values of baryon branching ratios into sterile neutrinos are propotional to 
squared values of corresponding active-to-sterile mixing angles. The latter are limited from 
above due to negative results of direct searches for sterile neutrinos. For mixing with electron 
neutrinos, the strongest limits are \6ex\'^ ^ 3 ■ 10~^ — 3 ■ 10"'' for rriK < Tn,^^ < 1 GeV from 
BEBC experiment [9], |^exp < 10"^ for 1.5 GeV < < 2 GeV from CHARM [lO] and 
l^exP < 1 ■ 10"^ - 1 • 10"^ for 2 GeV < < 3 GeV from HRS [llj. For mixing with muon 
neutrinos the limits are l^^^p < 5-10-^-l-10-^ for <m^^<l GeV, < e-lO-^-l-lO"^ 

for 1.5 GeV < m^^ < 2 GeV from NuTeV [12], l^^^p < 3- 10-^-4- 10"^ for 2 GeV < m^, < 2.5 
GeV from CHARM II experiment [I3] and \d^x\^ < 5 ■ 10"^ - 1 ■ 10"^ from HRS experiment [H]. 

Consequently, in models with sterile neutrino mass ranging within ttlk ^ rn^^ < 1 GeV one 
estimates from Fig. [1] and the above limits: 

Sr(Ae A + e+ + z/,) < 2.1 ■ 10-^ 
Br{A, ^ A + /i+ + z/,.) < 3.5 ■ 10'^ 

In models with 1.5 GeV < rriu^ < 2 GeV: 

Br{Ah ^ Ae + e- + u,) < 2.8 ■ 10"^ - 1.3 ■ 10"^ 
Br{Ah ^A, + i2- + z/,) < 1.7 ■ 10"^ - 1.3 • 10"^, 

where larger value on the right hand side corresponds to smaller value of and vice versa. 
In models with 2 GeV < m^^ < 2.5 GeV: 

Br{Ab ^ Ae + e- + z/,) < 1.3 • 10"^ - 1.7 ■ 10~^, 
Br{A, ^A, + fi- + V.,) < 3.9 ■ 10^^ - 1.5 • 10"^ 

Finally, in models with heavy sterile neutrinos in the range 2.5 GeV < m,y^ < 3 GeV: 

Br{Ah ^A, + e' + v,) < 1.7 ■ 10-^ 

and the same limit for the process A;, — A^ + ^~ + z/j.. 

Thus, to search for sterile neutrinos, one has to collect statistics of at least a few million 
heavy baryons. 
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